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Abstract. We describe a code designed for hot media (T > a few 10^ 
K), optically thick to Compton scattering. It computes the structure of 
a plane-parallel slab of gas in thermal and ionization equilibrium, illumi- 
nated on one side or on both sides by a given spectrum. This code has 
been presented in a previous paper (Dumont, Abrassart & Collin 2000), 
where several aspects were already discussed. So we focus here mainly on 
the recent developments. Presently the code solves the transfer of the con- 
tinuum with the Accelerated Lambda Iteration method (ALI) and that 
of the lines in a two stream Eddington approximation, without using the 
local escape probability formalism to approximate the line transfer. This 
transfer code is coupled with a Monte Carlo code which allows to take 
into account direct and inverse Compton diffusions, and to compute the 
spectrum emitted up to MeV energies, in any geometry. The influence 
of a few physical parameters is shown, and the importance of the den- 
sity and pressure distribution (constant density, pressure equilibrium, or 
hydrostatic equilibrium) is stressed. Recent improvements in the treat- 
ment of the atomic data are described, and foreseen developments are 
mentioned. 



1. Introduction 



The average continuum observed in Active Galactic Nuclei (AGN) is shown in 
Fig. 1 from the optical to X-ray range, as obtained from a composite optical- 
UV spectrum of Francis et al. (1991) and Zheng et al. (1997), and from an 
average (slightly modified) soft X-ray spectrum of Laor et al. (1997). It is 
widely admitted since several years that the optical-UV part of this continuum 
is produced by a warm (temperature of the order of 10^ — 10^ K), optically 
thick or effectively thick medium, most likely an accretion disc. Moreover, one 
deduces from the spectral distribution in the hard X-ray range, in particular 
the 30 keV hump observed in many Seyfert 1 galaxies, and from the presence of 
other features like the iron K line around 7 keV, that this medium is irradiated 
by a hard X-ray source, which is partly reprocessed and reemitted (improperly 
named "reflected"). So one sees a combination of the primary and of the re- 
flected spectra (Ross & Fabian 1993, and many subsequent works). In another 
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Figure 1. Average continuum of radio loud and radio quiet Active 
Galactic Nuclei, in arbitrary units. 



interpretation the UV-X continuum is produced by clouds from a disrupted disk 
(Collin-Souffrin et al. 1996, Czerny & Dumont 1998), and in this case one would 
also see the spectrum emitted by the non illuminated surface of the clouds. In 
both interpretations the UV-soft X spectrum is emitted by a dense and thick 
shell (density higher than 10^^ cm^^, Thomson thickness higher than unity), 
irradiated by an X-ray continuum. Similar conditions are met in X-ray binary 
stars. 

It is necessary to compute the structure of such an irradiated medium to 
predict the observed spectrum. Therefore one has to build "intermediate" codes, 
between those designed for planetary nebulae or Narrow or Broad Line Regions 
in AGN, and those designed for stellar atmospheres. These codes have to be 
valid for thick or semi- infinite dense media, eventually hot, irradiated by a non 
thermal continuum extending in the hard X-rays, to cover the whole range of 
situations encompassed in AGN and in binary stars. 

Owing to the high optical thickness of the medium in several frequency 
ranges, such codes require that the transfer of both the continuum and the lines 
be solved in an "exact" way, that is avoiding approximations such as local escape 
approximation (for the lines) or one stream approximation (for the continuum). 
Since the medium is generally dense and sometimes close to Local Thermody- 
namical Equilibrium (LTE), they require that all processes and inverse processes 
be carefully handled. Being irradiated by an X-ray continuum, the medium con- 
tains a large number of ionic species, from low to high ionization, which should 
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all be introduced in the computation. Finally, the medium being hot and thick, 
not only Thomson, but also Compton scattering, should be taken into account. 

We have undertaken to build a code which satisfies these requirements, to 
compute the emission spectrum produced by irradiated thick and hot media 
like those commonly assumed in AGN close to the black hole, in a wide range 
of photon energy. Precisely wc have built several interconnected codes, which 
allow more flexibility. The ensemble is far from being perfect and still contains 
several approximations which restrict its use, but we improve it gradually. 

TITAN is designed to study the structure of a warm or hot thick photoion- 
ized gas, and to compute its emission - reflection - transmission spectrum from 
the infrared up to about 20 keV. It solves the energy balance, the ionization and 
the statistical equilibria, the transfer equations, in a plane-parallel geometry, for 
the lines and the continuum. Then, given the thermal and ionization stratifica- 
tion, the computation of the emitted spectrum from 1 kcV to a few hundreds 
keV is performed with NOAR which uses a Monte-Carlo method taking into 
account direct and inverse Compton scattering (it allows also to study various 
geometries). In a previous paper (Dumont, Abrassart & Collin 2000) TITAN 
and its interconnection with NOAR were described. We recall here the main 
characteristics of the code, and we describe recent improvements as well as new 
results, focussing only on the aspects which are not treated in a standard way. 

We briefiy summarize below the physical processes (Section 2) , the transfer 
method and the iteration procedure (Section 3). The infiuence of the physical 
parameters, of the atomic data, and of the density distribution, are discussed in 
Section 4. The coupling of TITAN and NOAR is briefly described in Section 5. 



2. Generalities about TITAN 
2.1. Physical processes 

In TITAN, the physical state of the gas (temperature, ionic abundances and level 
populations of all ionic species), is computed at each depth, assuming stationary 
state, i.e.: local balance between ionizations and recombinations, local balance 
between excitations and deexcitations, local energy balance, and finally total 
energy balance (equality between inward and outward fluxes). 

The ionization equilibrium equations include radiative ionizations by con- 
tinuum and line photons, collisional ionizations and recombinations, radiative 
and dielectronic recombinations, charge transfer by H and He atoms, the Auger 
effects, and ionizations by high energy electrons arising from ionizations by X- 
ray photons. The emission-absorption mechanisms of the continuum include 
free-free and free-bound processes, two-photon process, and Thomson electron 
scattering. Inverse processes (except autoionization for dielectronic recombina- 
tions, and recombinations onto excited levels for some ions, see next section) 
are computed through the equations of detailed balance. Induced processes are 
taken into account. Energy balance equations include the same processes and 
Compton heating/cooling. 
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2.2. Atomic data 

Hydrogen and hydrogen-like ions are treated as 5 or 6-lcvcl atoms. In a first 
step, in order to save computation time, non hydrogen-like ions were treated 
with a rough approximation: interlocking between excited levels was neglected 
and the populations of excited levels were computed separately using a two-level 
approximation. This approximation does not predict correctly the details of the 
line spectrum, including the resonance lines themselves. Moreover ionizations 
from excited levels were not taken into account, though recombinations onto 
excited states were taken into account in the ionization equilibrium. This ap- 
proach is is not correct close to LTE, and even far from LTE, if ones wants to 
predict accurately detailed spectral features (cf. Section 4.2). However it gives 
a correct overall spectral distribution in this last case. Such arc the upper layers 
of an irradiated atmosphere where at the same time the fractional abundances of 
heavy elements in high ionization states are important, and the ground level of 
these ions is strongly ovcrpopulatcd with respect to LTE. But it is not the case 
in the deep layers of a disk atmosphere, emitting the UV part of the spectrum, 
where the density is high, and the influence of the external irradiation is small. 

This is why helium-like ions treated as complete 8-level atoms and lithium- 
like ions treated as complete 5-level atoms have been recently implemented in 
the code. All processes are thus taken into account for each level. This treatment 
is valid as far as the upper levels are in LTE with the continuum. To handle 
intermediate cases, we replace the missing upper levels which are very important 
especially for helium-like ions by adding an additional recombination rate, not 
compensated by ionization. 

In the future we plane to add several other levels, and to sum the contri- 
butions of the higher levels to handle at the same time cases close to and cases 
far from LTE. 

Photoionization cross sections are fitted from Topbase (cf. Cunto et al 
1993). Cross sections of neutral and once ionized heavy elements are not well 
represented, but they do not have any influence on these hot media. 

The gas composition include 10 elements (H, He, C, N, O, Ne, Mg, Si, S, 
Fe), and all their ionic species. 

The slab is divided in about 300 layers with variable thicknesses. The 
code can work with a prescription for the density or for the pressure, such as a 
constant pressure or a pressure corresponding to hydrostatic equilibrium. 

2.3. Illumination 

The external illumination can be characterized by an "ionization parameter", 
which we choose here as ^ = 47rFinc/?^Hj -^inc being the frequency integrated 
flux incident on one side of the slab and tt-h the hydrogen number density at this 
surface. Note that there are several other definitions of the ionization parameter, 
for instance U, or S, but they correspond basically to the same concept, namely 
a parameter which determines the ionization state. It is possible to take also 
into account a fiux incident on the back side of the slab. 

Presently we mimic a semi-isotropic illumination, but with the ALI transfer 
method which is implemented (cf. below) it will be possible to take into account 
the angular distribution of the incident radiation. 
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Figure 2. Reflected spectrum, with (full line), and without (dashed 
lines) ALI transfer for the continuum, for the reference model defined 
in Section 4; dotted line: incident spectrum; the spectrum is displayed 
with a spectral resolution of 30. 

3. Radiation transfer 

We are concerned by media with a continuum optical thickness larger than unity 
in a range of frequencies rich in emission or in absorption lines, and having a 
very inhomogeneous structure. Moreover in the hot photoionized layers, the 
absorption coefficient is weak at all wavelengths and Thomson diffusion is dom- 
inant. On the contrary absorption dominates and must be carefully treated to 
compute the spectrum emitted by the non illuminated cold part of the medium. 
Radiation transfer requires therefore some attention. 
The transfer equation is 



where z is the distance to the illuminated surface, /i is the cosine of the angle 
between the normal and the light ray, k is the absorption coefficient (for the 
continuum it is due to photoionizations and free- free transitions), a is the diffu- 
sion coefficient (here Thomson scattering), r] is the emissivity (due to radiative 
recombinations and free- free emission), and J is the mean intensity, equal to 
/ (all these quantities are local, / depends on the direction and on the 

frequency u, and J, k and r] depend on the frequency). 

Introducing the source function S and the optical depth r, Eq. 1 can be 
written: 




(1) 
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,f = -I + S (2) 



with 

aJ + 1] 



(3) 



a + K 

which can be separated between J-dependent and J-independent components: 

S=il-e)J + eB. (4) 

The formal solution is: 

Ar) = \jdlij S{t)e-^'-^y^dt/ii. (5) 
3.1. Transfer of the continuum: the ALI method 

In the previous version of the code (Dumont et al. 2000) wc used the same 
method for line and continuum transfer, described in Section 3.2. Presently 
the transfer of the continuum is performed with the Accelerated Lambda 
Iteration method (ALI), which was set up with the help of subroutines kindly 
provided to us by F. Paletou (cf. Auer &: Paletou 1994, and Kunasz & Auer 
1988). 

In this method J(t) is written with the help of an operator A: 

J{r) = K[S{t)]. (6) 

ALI uses an operator splitting, described first by Cannon (1973), then by Olson 
et al. (1986): 

A = A* + (A-A*) and S = S* + 5S (7) 

where A* is an approximate A operator chosen to facilitate the inversion, and 
S* has been calculated in the previous iteration. One obtains: 

(5S = [l-(l-e)A*]-i[5-5*]. (8) 

We use here a parabolic interpolation in the evaluation of the integral giving 
the intensity from the source function, and an acceleration type Ng (in this 
method, at every fourth iteration a least squares extrapolation is made in order 
to reduce the residuals, cf. Ng 1974). 

To perform this computation, one needs to know the opacities and emis- 
sivities at each frequency as functions of z, hence the temperature and the 
populations at each depth, so an iteration procedure is required. 

• For each layer, starting from the illuminated side, the ionization and ther- 
mal balance equations are solved by iteration; the temperature, the opacity 
and the emissivity are computed; 

• when the back side of the cloud is reached, the transfer is solved with the 
ALI method, and Ju{z) is computed; 
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Figure 3. Profiles of I^, and for OVIII Lya, for the reference 
model defined in Section 4. The curves are labelled by the value of 
z/H. 
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• the whole calculation is repeated until convergence with the new values 
of Jv{z). It is stopped when the energy balance is achieved for the whole 
slab (i.e. when the flux entering on both sides of the slab is equal to the 
flux coming out from both sides). This method is almost as slow as the 
previous one, albeit more secure. 

Fig. 2 shows that the difference in the spectrum computed with the previous 
transfer method and with ALI is very small. It is certainly due to the fact that 
the continuum optical thickness is never very large (at most 10^). This is not 
the case at the center of intense lines, which are typically 100 times thicker. We 
are presently implementing ALI for the line transfer and we expect to find more 
important differences for a few lines like OVIII Lq which are not completely 
converged even after 1000 iterations (cf. Dumont et al. 2000). 

3.2. Line transfer 

Since ALI is presently not implemented for the line transfer, it is treated like in 
the previous version with a simple method based on the Eddington two-stream 
approximation. In this approximation the transfer equations are written: 



TT = -{l^^+nK+nlu+^u (9) 



V3 dz ' " ' 2' " '2' 



V3 dz ' " 2' 2' 

The mean intensity Ji, is equal to (/j" + /j7)/2, and the flux F,^ defined as 
usual by / ly cos 6 duj is equal to (/+ — I^) 27r/\/3. The "reflected" fiux is equal 
to I^{0) 27r/\/3, and the outward flux to I^{H) 27r/\/3. The optical depth and 
the total optical depth are defined as: 

r^(z) = rVS{K^ + a)dz' and % = t^{H) (10) 
Jo 

where H is the total thickness of the slab. For I^{z) the formal solution of the 
transfer equations between z — dz and z is: 

I+(z) = I+{z - 5z)e-^^^ + e"^- S^{t)e+^dt. (11) 

We approximate the source function Si,{t) by a constant in the interval 
z,z + dz. It leads to: 



1 - e 



I+{z) = I+iz - 5z)e-'-" + [S,{z - 5z) + S,{z)] (12) 

and: 

I-{z) = I-{z + 5z)e-^^'^ + ^—^^[S^^iz + 5z) + S^{z)]. (13) 
We calculate Si,{z) and the mean intensity Ji,{z) from Eq. 3 and from: 
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Figure 4. Temperature versus depth for different values of ^ in erg 
cm s~^, the other parameters being the same as in the reference model. 



assuming a constant 5z (actually it is not constant so the formulae are more 
complicated but do not deserve to be given here). 
The procedure is then the following: 

• for each layer, after computation of the temperature and the source func- 
tion, I^{z) is transferred through the layer according to Eq. 12, while 
Sy[z) and Ty{z) are buffered for each layer and each frequency; 

• when the back side of the cloud is reached, new values of I~ {z) are cal- 
culated, starting from the back side where I~{H) is given by Eq. 12, and 
using Eq. 13 with the previous values of Si,{z) and Ty{z) (actually we use 
the values given by several previous iterations to accelerate the procedure 
with a Ng method). 

This method requires a large number of iterations to get a complete conver- 
gence of the model, and one cannot be sure to have reached the correct value. 
It is why we are now changing for the ALI method. 

We assume presently complete redistribution in the lines, and partial redis- 
tribution in intense lines is mimicked by assuming complete redistribution in a 



Ju{z) 



I+{z-5z)+I-{z + 5z) . 



(14) 
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Figure 5. Reflected and outward spectra for different values of ^ in 
erg cm s~^, the other parameters being the same as in the reference 
model; the spectra are displayed with a spectral resolution of 100. 

Figure 6. Reflected and outward spectra for different values of the 
column density in cm~^, for ^ = 100 erg cm s~^, the other parameters 
being the same as in the reference model; the spectra are displayed 
with a resolution of 100. Dotted line: the incident spectrum. 

pure Doppler proflle. It will be possible to treat partial redistribution with the 
ALT method. 

A line proflle is represented by a Voigt function at several frequencies dis- 
tributed symmetrically around the line center. The frequencies in the proflle are 
chosen to take into account the fact that the Doppler width may decrease by one 
order of magnitude inside the slab, owing to large temperature variations. More- 
over, to solve the statistical equilibrium equations one needs to know the total 
line intensity weighted by the proflle / Ji,4>{i')dv, while to compute the thermal 
and ionization equilibria one simply needs the line intensity integrated over fre- 
quencies / Judv- For intense lines, the flrst integral is dominated by the Doppler 
core, and the second by the wings. The set of frequencies chosen to describe the 
line must therefore cover at least two orders of magnitude. The integrals over 
the line profiles are achieved with a 15-point Gauss-Legendre quadrature. We 
also assume that the lines do not overlap, which is valid as far as multiplets are 
treated globally, and because the number of lines taken into account in TITAN 
is relatively small (presently 580, including subordinate lines in the optical-UV 
range) . 

As shown by Fig. 3, the profiles of /~ and for OVIII Lya vary consid- 
erably from the surface to the deepest layers. 

4. Some results of TITAN 

We shall call "reference model" a model with the following characteristics: 

• it is a plan-parallel slab with a constant density, equal to 10^^ cm~^, 

• its column density equals 10^^ cm~^ (corresponding to a Thomson thick- 
ness Tes of 80), 

• it is illuminated on one side by an incident power law continuum propor- 
tional to from 0.1 eV to 100 keV, 

• the ionization parameter is ^ = 10^ erg cm s~^, 

• there is no illumination on the other side, 

• the abundances are (in number): H: 1, He: 0.085, C: 3.3 10"^, N: 9.1 10"^, 
O: 6.6 10-^, Ne: 8.3 10-^ Mg: 2.6 10-^ Si: 3.3 10-^ S: 1.6 10-^ Fe: 3.2 
10-^ 
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In several examples given below, we also use a smaller value of the ionization 
parameter, ^ = 10^ (corresponding to a value of the usual ionization parameter 
U equal to 0.88), allowing an easier comparison with computations made with 
CLOUDY or with XSTAR. 

Figs. 4 and 5 show the results corresponding to different values of the ioniza- 
tion parameter, the other parameters being the same as in the reference model. 
These results were obtained with the previous version of the code. Indeed each 
model takes a long time to run, and all models have still not been recalculated 
with the ALI method. However, we have checked in several cases that the results 
obtained with the new version differ very slightly from the previous ones. 

Fig. 4 displays the temperature versus the depth. Close to the surface the 
temperature is high and reaches the Compton temperature for the highest value 
of ^. The extension of this region (called the "hot skin") increases with the value 
of ^. Its optical thickness reaches Tcs ~ 10 for ^ = 10^. As we shall sec later, 
this is quite different from the case of constant pressure, where the thickness of 
this zone is always limited to Tes < 1. This hot skin is almost absent for low 
values of ^. Below the hot skin the temperature decreases smoothly, again at 
contradistinction with the case of a constant pressure. This decrease is followed 
by a region of a quasi constant temperature, which is in a state of "quasi-LTE" , 
i.e. each process important for the energy equilibrium (lines, free-bound from 
and onto ground levels, free-free) is exactly balanced, due to the high optical 
thickness and relatively high density. However the region is not in LTE, in the 
sense that the fractional ionic abundances do not satisfy Saha equations. We 
shall discuss this issue in more detail in the next section. 

On Fig. 5 are displayed the "reflected" spectrum emitted by the illuminated 
side of the slab, and the "outward" spectrum emitted by the back side. Note 
that the outward spectrum is distinct from the transmitted one, as it includes 
also emission (actually our transfer method does not allow to separate emission 
and transmission, contrary to approximate methods using for instance "on the 
spot" approximation). The outward spectrum lies only in the optical near UV 
range, owing to the low temperature of the back side (cf. Fig. 4) . This is due to 
the large column density of the model which does not allow any transmission in 
the X-ray range. Note the strong edges in emission (including a Balmcr edge), 
and the cut off of the spectrum at a few tens of eV, due to absorption by Hell 
and by low ionization heavy element species. The reflected spectrum shows a 
characteristic shape with a smooth "powerlaw" soft X-ray excess, and displays 
many lines and ionization edges whose intensity decreases with increasing values 
of ^, owing to the increasing influence of pure reflection when the ionization 
state increases. Note that the Lyman discontinuity is in emission for low ^ and 
in absorption for high ^. 

Fig. 6 illustrates the influence of the column density on the reflected and 
on the outward spectra for ^ = 100, the other parameters being the same as in 
the reference model. First we see that the reflected spectrum "saturates" above 
10^^ cm~^, because all the available photons have been absorbed and no longer 
succeed in escaping from the illuminated side. Second, the optical-UV part of 
the outward spectrum saturates above 10^^ cm~^, but not the X-ray part, which 
is due to transmission. As the column density decreases, the slab becomes more 
transparent and less emissive, so the outward spectrum becomes basically only 
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the transmitted spectrum, while the reflected spectrum decreases in intensity. 
Note that the "saturation" of the reflected spectrum would occur at a larger 
value of the column density for a larger ^, as the hot reflecting skin becomes 
thicker. It is thus necessary to take into account deeper layers, exceeding a 
Thomson thickness Tcs of few units, to compute correctly the reflected spectrum 
in the UV range for £, > 10^, when dealing with constant density slabs. This has 
always been overlooked in the past. 

4.1. Influence of the excited levels 

In the first version of the code, non hydrogen-like ions were treated with a rough 
approximation: the populations of the excited levels were computed separately 
using a two-level approximation, i.e. interlocking between excited levels was 
neglected, as well as ionizations from these levels. Moreover only a very small 
number of levels (at most four) were taken into account. However recombinations 
onto excited states were taken into account in the ionization equilibrium, and 
the transfer of these photons was treated in an approximative way, according to 
an escape probability formula proposed by Canfield & Ricchiazzi (1980). Each 
recombination was also assumed to produce a resonant photon after cascades. 

Hydrogen and hydrogen-like ions were - and are - treated as 6- level atoms. 
Levels 2s and 2p are treated separately, while full 1-mixing is assumed for higher 
levels. All processes including collisional and radiative ionizations and recombi- 
nations are taken into account for each level (cf. Mihalas 1978). In the previous 
version, recombinations onto levels higher than five were not taken into account, 
which amounted assuming that the higher levels were in LTE with the contin- 
uum. It was true for hydrogen in the conditions in which this code is presently 
used, but not for other hydrogen-like ions. We have therefore recently added 
recombinations onto higher levels, and we intend to implement "fictive" levels 
to take ionizations from upper levels into account as it is done for instance in 
CLOUDY or XSTAR. Level populations are obtained as usual by matrix inver- 
sion. 

The rough treatment of non hydrogen-like ions could have several conse- 
quences. 

• First it should not predict correctly the details of the line spectrum, since 
it neglects subordinate lines and even important resonance lines. 

• Also it can lead to a bad estimation of the energy losses and consequently to 
a change of the temperature. The recombination continuum corresponds 
generally to a smaller absorption coefficient than the neglected lines, and 
therefore can escape more easily from the deep layers. It would corre- 
spond to an overestimation of the losses and an underestimation of the 
temperature. But other parameters, such as the ionization state, play also 
a crucial role in the energy balance, so the net effect is difficult to predict. 
Fortunately in these models cooling due to lines is not predominant except 
in some regions. 

• It can lead to important errors in the fractional ion abundances. This effect 
is "self-cumulating", as an increased proportion of a given ion leads to a 
stronger ionization edge, which in turns inhibits further photoionizations. 
So again the net effect is not predictable. 



APS Conf. Ser. 



13 




Figure 7. Comparison of the temperature computed with (full line) 
and without (dashed lines) full eight-level He-like and full five-level Li- 
like ions, for the reference model, versus the distance to the illuminated 
side. 
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• Finally recombinations onto excited levels are not balanced by ionizations 
from these levels. It is not important in conditions prevailing in a pho- 
toionized slab far from equilibrium (like for instance in the hot skin of 
an irradiated accretion disk), as the lowest atom levels are strongly over- 
populated with respect to LTE. But this treatment is not valid in the 
innermost layers of a disk atmosphere, where the density is very high so 
complete LTE is almost reached. 

We have recently introduced interlocked multi-level computations for He- 
like and Li-like ions, in order to obtain a better description of the spectrum, and 
to treat cases closer to LTE. However we are for the moment restricted to 8 levels 
for He-like ions and 5 for Li-like ones, which is not sufficient close to LTE. We 
take into account recombinations onto the levels higher than 8 for Hc-likc ions 
(respectively 5 for Li-likc), adding them as a contribution to the recombination 
of the highest level. This treatment is similar to the old one, except that it 
takes into account many subordinate lines and ionizations from excited levels, 
and it is valid as far as LTE is not reached for the highest levels and ionizations 
are dominated by the ground level. Like for hydrogen-like ions we plane to add 
several other levels, and to mimic the contributions of the higher levels. 

We mention also that radiative recombinations are computed differently 
according to the relative values of kT to the photon energy bin, in order to get 
an accurate frequency dependence. 

Figs. 7 to 10 illustrate the influence of the new treatment with respect to 
the old one, for the reference model (except that the thickness is only Tes = 25 
instead of 80; we have already mentioned that above this value nothing varies 
anymore for ^ = 10^, cf. Fig. 4). In this model the level populations and 
the ionization degrees are far from LTE, even in the deepest layers. It is quite 
comparable to a classical photoionized model, where excited levels are populated 
by radiative recombinations and depopulated by deexcitations onto lower levels. 
Note that these results are obtained with the new transfer treatment. 

Fig. 7 shows the new temperature profile. Except at the surface T is 
increased with respect to the previous computation. The change is quite impor- 
tant, reaching a factor two in the transition layer. Fig. 8 shows the change of 
fractional abundances of Hydrogen, Helium, and Iron. It is relatively important 
in the hot skin for Hel, but the abundance of this ion is completely negligible in 
this region, so it should not induce any change in the emitted spectrum around 
the ionization edge, i.e. 20 eV. The result of a computation with full He-like 
and Li-like atoms, but without additional recombinations on the highest level, 
is also shown. Except precisely for Hel in the hot skin, the results with and 
without the additional term are quite similar. 

Figs. 9 and 10 display the reflected spectrum for this model, computed with 
the old and the new atomic treatment, with a low resolution in the whole energy 
range (Fig. 9), and with a higher resolution in a smallest but interesting energy 
range (Fig. 10). 

The spectral distribution docs not differ much for the two treatments. This 
is not the case of the detailed features, as it can be seen on Fig. 10. The flux of 
the main lines of H-like ions in the reflected spectrum are almost not changed 
(which is expected since the treatment is the same). For the He-like and Li-likc 
ions, the total flux for all lines is multiplied by a factor 3 or 4. To illustrate this 
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Figure 8. Fractional abundances of Hydrogen, Helium and Iron, com- 
puted with (dotted line) and without (full lines) full eight-level He-like 
and full five-level Li-like ions, for the reference model, versus the dis- 
tance to the illuminated side. The result of a computation with full 
He-like and Li-like atoms, but without additional recombinations on 
the highest level, is also shown (dashed lines). 
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Figure 9. Reflected spectrum computed with and without fuh eight- 
level He-like and full five-level Li-like ions, for the reference model. The 
result of a computation with full He-like and Li-like atoms, but without 
additional recombinations on the highest level, is also shown. The 
spectra are vertically shifted for a better reading. They are displayed 
with a resolution of 30. A few intense lines are identified. 
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Figure 10. Reflected spectrum computed with and without full eight- 
level He-like and full five-level Li-like ions, for the reference model. The 
result of a computation with full He-like and Li-like atoms, but without 
additional recombinations on the highest level, is also shown. The 
spectra are vertically shifted for a better reading. They are displayed 
with a resolution of 100. A few intense lines are identified. 
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effect, Fig. 11 displays the fluxes of the OVI, OVII, and OVIII lines. Even the 
resonant lines differ in both treatments, as the energy is distributed among a 
much larger number of lines: for instance the OVIII line at 665 eV decreases by 
one order of magnitude with the new treatment. The net increase of line flux is 
simply due to the redistribution of continuum recombination photons into line 
photons through cascades, an effect which was only partly taken into account in 
the previous treatment. 

Besides, we have added the result of a computation performed with the 
complete treatment of He-Li-like ions, when suppressing the recombination rates 
added on the highest level (cf. Figs. 9 and 10). This suppression does not affect 
at all the spectral distribution, and only slightly the detailed features of the 
spectrum. 

The next further step is to introduce forbidden lines, which could be im- 
portant for the cooling of the medium in some places. 

Note that all the results given in the following sections are obtained without 
Li-like and He-like complete atom models. 

4.2. Influence of the density distribution 

Constant pressure : 

Up to now, we have presented results only for slabs of constant density. It is 
a peculiar the pressure is thus varying according to the density and the 

temperature. When the pressure is imposed and not the density, one is faced to 
the well-known problem of thermal instability described first in the framework 
of AGN by Krolik, McKee & Tarter (1981). If the pressure is imposed, three 
(or even more) solutions to the Gain=Loss equation can exist, and one of these 
solutions is unstable. The gas has therefore to shift either to the Compton 
dominated hot solution at a temperature close to the Compton value, or to the 
cold solution dominated by atomic processes, at a few 10^ K. 

Krolik et al. were considering optically thin media with a uniform temper- 
ature, whose value is given by the external irradiating spectrum. Here we are 
considering a stationary case, where the temperature changes across the slab 
because the mean intensity Jj^ is modified as a function of depth. Since is 
absorbed in the soft X-ray band, its spectral distribution becomes more dom- 
inated by UV as the depth increases, and such a spectrum docs not lead to 
the S-shape curve characteristic of thermal instability (cf. Krolik et al. 1981). 
However, even if there is no instability, the temperature decreases very rapidly 
with increasing depth, as the constancy of the pressure leads to an increase of 
the density, therefore to an increase of the cooling, and to a further decrease of 
temperature (note that the temperature is not strictly inversely proportional to 
the density, as radiation pressure generally exceeds gas pressure). 

This is illustrated on Fig. 12 which shows the Gain-Loss function versus the 
temperature for different layers. The model is a slab with a constant pressure, 
a density at the surface equal to 10^^ cm~^, an ionization parameter equal to 
1000. The illuminating spectrum is a power law F^, oc u~^''^, chosen because such 
a flat spectrum is well known to correspond to an "S-shape" curve with a strong 
instability. However we see that as the temperature decreases with the depth, 
there are never multiple solutions of the Gain=Loss equation. Note that it is by 
no means a general result. In particular, a medium with a more rapid variation of 
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Figure 11. Line fluxes with (open symbols) and without (full sym- 
bols) full eight-level He-like and full five-level Li-like ions, for the ref- 
erence model. OVL circles; OVIL losanges; OVIIL squares. 
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Figure 12. Gain-Loss function versus the temperature, for different 
layers across the slab, for a constant pressure model. 



the pressure with the depth, like given by hydrostatic equilibrium, could likely 
lead to a thermal instability, as the transition to the cold region occurs in a 
region where the spectrum is less absorbed than in a constant density medium 
(cf. Nayakshin, Kazanas & Kallman 2000). Also if the UV band is absorbed 
before the soft X-ray one, which is the case when the ionization parameter is 
small, it should lead to a thermal instability. 

Fig. 13 displays the temperature versus depth for slabs of constant density 
and slabs of constant pressure, with different ionization parameters. Note that 
it would be better to speak of "ionizing flux" instead of ionizing parameter in 
the constant pressure case. Actually the comparison is performed with the same 
incident fluxes. As expected, and in contrast with the constant density case 
where there is a smooth decrease of the temperature with increasing depth, the 
constant pressure case gives a sharp transition between the hot and dilute, and 
the cold and dense, parts of the slab. Moreover the thickness of the "hot skin" 
is smaller for the same flux in the constant pressure case, in particular for high 
illumination. For instance, for ^ = 10^ the Thomson thickness of the hot skin is 
about unity in the constant density case, and only 0.4 in the constant pressure 
case. All these effects have strong consequences on the emission spectrum. 

For constant pressure the spectrum displays a narrow UV bump, while the 
constant density the UV bump is more extended in frequency, as it corresponds 
to intermediate temperatures which do not exist in the constant pressure case 
(cf. Fig. 14). The difference between the two cases is particularly apparent for 
high illumination, as expected. It is interesting to note that the Lyman edge 
is quite weak in the constant pressure case, since the UV spectrum is closer to 
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Figure 13. Temperature versus depth for different values of in erg 
cm s~^, for constant pressure and constant density. The constant den- 
sity slabs have a density of lO^'^ cm~^, and the density at the surface of 
the slab with constant pressure is equal to the same value. The other 
parameters are the same as in the reference model. 
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Figure 14. Reflected spectra computed with constant pressure (full 
line) and with constant density (dashed lines) for the models of the 
previous flgure; the spectra are displayed with a resolution of 30. 
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a blackbody. This is a point that constant density models have difficulties to 
account for in AGN. 

The distribution of fractional ionic abundances is also completely modified. 
As an illustration, Fig. 15 displays the fractional abundances of iron versus 
depth. The most remarkable difference is the presence at t^s ~ 1 of low ionization 
species, such as FcIV or FeV, which are completely absent in the constant density 
case. One notes also that the extension of highly ionized species is smaller than 
in the constant density expected. For instance FcXIV is present up 

to Tes — 0.4 for constant pressure, instead of 1.2 for constant density. As a 
consequence, the Iron K lines produced are different from the constant density, 
as we shall see later. 

Hydrostatic equilibrium : 

In order to describe an irradiated accretion disk, we have introduced in 
collaboration with A. Rosanska and B. Czerny a variation of the density given 
by the hydrostatic equilibrium law. TITAN treats then only the layer at the 
surface of the disc where the diffusion approximation does not hold. In this case 
there is a flux incident on the back side of the slab, equal to the flux generated 
by viscous dissipation inside the disk. Iterations are needed with the hydrostatic 
equilibrium equations (for more details see Rozahska et al. 2001). 

Fig. 16 shows the results for the following case: Mbh = 10^ M©, m = 0.03, 
R = 10-Rsch) and an X-ray irradiating powerlaw spectrum, with a spectral index 
0.9 in flux, and an X-ray flux Fx = 3.5 x 10^^, i.e. equal to half the viscous 
flux. Note that the temperature decreases abruptly like in the constant pressure 
case, and that the hot skin is very thin. The temperature reaches that given 
by viscous dissipation at about Tes = 0.4. The reflected spectrum displays a 
prominent blue bump with no soft X-ray excess, and a strong absorption in the 
UV range due to low ionization species in the deep layers of the atmosphere. 
It is important to note that such a spectrum differs strongly from that of an 
irradiated disc in which the density is assumed to be constant (Ross & Fabian 
1993, and subsequent works, cf. Rozahska et al. 2001 for a detailed discussion). 
Our results are similar to those of Nayakshin, Kazanas & Kallman (2000). 

5. Coupling of TITAN and NOAR 
5.1. Monte Carlo code NOAR 

In order to compute the effect of Compton scattering on the high energy part 
of the reprocessed radiation and to determine precisely the Compton heating- 
cooling rate in the energy balance equation, one needs another numerical ap- 
proach. In this aim Abrassart (1998) developed a Monte Carlo code, NOAR 
(cf. Dumont et al. 2000). The asset of such an approach is also that it enables 
to investigate an arbitrary geometry and to determine the angular dependence 
of the observed spectra. Moreover, it allows to easily extract time variability 
information. 

NOAR takes into account direct and inverse Compton scattering, according 

to the method proposed by Pozdniakov, Sobol & Sunyacv (1983) and by Gorecki 
&; Wilczewski (1984). Given all the fractional abundances and the temperature 
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Figure 15. Fractional abundances of iron versus depth, in the case of 
constant pressure or constant density slabs, for the incident flux 8.10^^ 
erg cm~-^ s~^ for one of the model of the previous figure = 1000). 
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Figure 16. Temperature versus depth, and reflected spectrum, ( 
tained for an illuminated accretion disk in hydrostatic equilibrium, 
described in the text. From Rozahska et al., 2001. 
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in each layer provided by TITAN, it computes the absorption cross sections in 
each layer. Free-free absorption is taken into account, as well as recombination 
continua of hydrogen and helium like ions (the spectral distribution of the re- 
combination continuum is determined by the local temperature). The proper 
yields to account for the competition with the Auger effect are used for the 
fluorescence lines. Fluorescence of Iron XVII to XXIII is suppressed by reso- 
nant trapping (but not in TITAN). For the pseudo fluorescence of Iron XXV 
and XXVI, a simplifying prescription is adopted which includes the two most 
probable outcomes of a K shell photoionization, that is: direct recombination 
on ground level or L-K transition. Other line emission is not included. 

One use of NOAR is to provide TITAN with the local Compton gains and 
losses in each layer. This is necessary, because Compton heating-cooling is 
dominated by energy losses of photons > 20 keV for high values of ^. In Dumont 
et al. (2000) it was shown that given an incident spectrum the dependence of 
the Compton heating-cooling over ^ ratio on the Thomson depth depends only 
slightly on ^ except for high values of ^. The Compton heating-cooling rate 
obtained with NOAR is fitted analytically as a function of z, which is transferred 
to TITAN. The spectra computed with TITAN are then completed, (cf. Fig. 
17, and also Figs. 14 and 16). 

5.2. Some results of TITAN and NOAR 

TITAN and NOAR lead to similar energy distributions of the reflected spectra in 
the 1-20 keV range. The reflected spectrum above 10 keV exhibits a "Compton 
hump" depending on the high energy cut-off of the incident spectrum, which 
is here 100 keV, and on the ionization parameter. The spectrum in the higher 
energy range is used for the energy balance. 

Comptonisation leads also to line broadening, and when it is important, to 
smearing of the ionization edges. Direct Compton downscatterings produce a red 
wing, while inverse Compton upscatterings produce a blue wing. The effect can 
be comparable to relativistic broadening near the black hole, and can account 
partly for the broad proflles of the iron line in AGN (Abrassart & Dumont 1998). 

Fig. 18 displays (with a relatively high resolution, in order to show detailed 
features) the spectrum around 7 keV in the constant density case, for different 
values of ^, the other parameters being as in the reference model. The feature 
near 7 keV is made of a mixture of iron edges and of several iron lines: 

- for ^ equal to a few hundreds erg cm s~^, the iron line is dominated by 
low ionization stages, smaller than iron XVII; the absorption edge is proeminent 
above 7 keV; 

- for ^ equal to 10^ erg cm s~^, FeXXV and FeXXVI lines are apparent, mod- 
ified by Compton scattering. The lines are significantly broadened, the broad- 
ening is asymmetric, the profile is skewed towards the red. The absorption edge 
is still important, but extends from 8 to 10 keV, due both to Comptonization 
and to the influence of several edges; 

- for ^ equal to 10^ erg cm s~^, the line has a large red Compton wing and 
a weak blue wing, and it is weak. The absorption edge is completely erased. 

Fig. 19 displays the iron line complex for the constant pressure case, to 
be compared to the constant density case. Both a "cold line" at 6.4 keV and 
FeXXV and XXVI lines are present in the constant pressure case, while only 
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Figure 17. Reflected spectrum for different values of ^ in erg cm s~^, 
the other parameters being the same as in the reference model; full 
lines: results of NOAR; dashed lines: results of TITAN. The spectra 
are displayed with a resolution of 100. 
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Figure 18. Reflected spectrum around 7 keV for different values of 
in erg cm s~^, the other parameters being the same as in the reference 
model; The broadening by Comptonization is important for ^ > 1000 
erg cm s~^. The spectra are displayed with a resolution of 200. 

highly ionized species are visible in the constant density case. For constant 
pressure the cold layers are indeed closer from the surface than in the constant 
density case, so the photons emitted by these layers are able to leave the medium. 



6. Foreseen improvements 

The coupling of TITAN and NOAR allows to compute the structure and the 
emission of hot Compton thick irradiated media, by solving consistently the 
energy balance of the medium. Moreover it allows also to take into account 
Compton broadening of the lines and of the edges in the X-ray range. 

In the previous description of the code (Dumont et al. 2000), we had already 
recalled the importance of the returning flux, which is often neglected, even for 
relatively low column densities, and we had compared the line transfer treatment 
for the lines with the escape probability approximation generally used in these 
problems. Here we have in addition stressed the influence of the pressure law, by 
showing the strong differences existing both in the structure and in the emitted 
spectrum, between a constant density medium, a constant pressure medium, 
or a medium in hydrostatic equilibrium. We have also shown that the atom 
models are of fundamental importance as they determine the ionization state, 
the energy balance, the temperature, and in fine the detailed spectral features 
which will be soon accessible through X-ray spectroscopy. 
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Figure 19. Reflected spectrum in the X-range in the case of constant 
pressure or constant density slabs, for the incident flux 8.10^^ erg cm~^ 
s~^. The spectra are displayed with a resolution of 200. 



All along this paper we have mentioned what improvements we intend to 
bring to the code. 

The Accelerated Lambda Iteration method is used now for the continuum. 
We will soon implement it for the lines, which are the most difficult to converge. 
It will in particular allow to treat partial redistribution in the lines. 

In collaboration with S. Coupe and M.C. Artru (Coupe et al. 2001) we have 
begun to take into account subordinate lines for all ions by solving a complete 
multi-level atom for Li-like and He-like ions, after H-like ions which were already 
introduced in the previous version of the code. In the future we plan to imple- 
ment other iso-electronic sequences, and to add to the present ones several other 
levels. In particular we will try to get a better representation of the highest 
levels for computations close to LTE. We will take into account forbidden lines, 
which could be important for the cooling of the medium in some places. We 
also intend to add to the line transfer an option allowing the UV and soft X-ray 
photons to escape from the medium through Compton scattering. 

Though there are still many improvements to perform, we presently use 
these codes to model AGN spectra in the UV and X-ray range, in conditions 
where they are valid. 
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